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Abstract. In Arsenic a pressure-induced solid-solid phase transition from the A7 into 
the simple cubic structure has been experimentally demonstrated [Beister et al. 1990 
Phys. Rev. B 41 5535]. In this paper wc present calculations, which predict that this 
phase transition can also be induced by an ultrashort laser pulse in As under pressure. 
In addition, calculations for the pressure-induced phase transition are presented. Using 
density functional theory in the generalized gradient approximation, we found that the 
pressure-induced phase transition takes place at 26.3 GPa and is accompanied by a 
volume change AV^ = 0.5 Og/atom. The laser-induced phase transition is predicted for 
an applied pressure of 23.8 GPa and an absorbed laser energy of 2.8 mRy/atom. 
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1. Introduction 

Femtosecond lasers create extreme nonequilibrium conditions in matter, namely, one 
with hot electrons and cold ions [1]. This is primarily because laser light interacts 
very strongly with electrons but not with ions, and secondarily because the electrons 
thermalize relatively slowly with the ions compared with the time needed for ultrafast 
structural changes (typically several ps vs several lOO's of fs). The possibility of heating 
the electrons in a material with an ultrafast laser pulse without immediately changing 
the temperature of the ions has paved the way to the discovery of interesting physical 
phenomena, such as, perhaps most importantly, laser-induced phase transitions that 
cannot be explained by the heating of the ions, but are caused by a change of the 
electronic bonding properties. 

In Arsenic a series of pressure-induced phase transitions exists. At ambient 
conditions As crystallizes in the A7 structure. At 25 GPa there is a transition to the 
simple cubic (sc) phase [2], followed by transitions [3] to the so-called As(III) structure at 
48 GPa and to body-centered cubic As at 97 GPa. The fact that As in the A7 structure is 
electronically stabilized by the Peierls mechanism |1] makes it a candidate for an ultrafast 
laser-induced solid-solid phase transition. This is of fundamental interest, because the 
majority of the experimentally studied laser-induced phase transitions shows ultrafast 
melting, for example, in Si |5], GaAs [6], Te [7], and InSb |8j. Examples of laser- 
induced solid-solid phase transitions are a disorder-to-disorder transition in amorphous 
GeSb [9], a metallic-to-semiconductor phase transition in SmS |TU], a monoclinic-to- 
rutile transition in VO2 [HI [l2], a ferromagnetic-to-paramagnetic transition in CoPta 
alloy films [13] , and an antiferromagnetic-to- ferromagnetic phase transition in FeRh [H] . 
The main topic of this paper is to analyze the possibility to induce a solid-solid phase 
transition by an ultrafast laser pulse in As under pressure. 

We now describe the relationship between the A7 and sc structures, the two phases 
of As between which we will predict a laser-induced transition (see section [s]). Starting 
from an sc atomic packing, the A7 structure can be derived in two steps, which are 
illustrated in figure [T} First the sc lattice is deformed by elongating it along one of 
the body diagonals (indicated by in figure The magnitude of this deformation 
is usually expressed by c/a, where c is the length of the lattice vector as and a is 
the length of ai (figure [l]). In the sc lattice c/a = ^/Q ~ 2.45. In the second step a 
Peierls instability causes the atoms to be displaced along the same diagonal, in opposite 
directions (figure [T|. The magnitude of the displacements is expressed by the atomic 
coordinate z: A value of 0.25 indicates no Peierls distortion. Deviations from z = 0.25 
give the magnitude of the displacement of the atoms in the as direction. In As at 
ambient pressure z = 0.228 ^15j, which leads to a displacement of the atoms amounting 
to 13% of the average distance between the planes [see figure [T|c)]. 
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Figure 1. Relation between the A7 and sc structures, a) Tlie sc structure. Tlie 
lattice vectors ai, a2, and belong to the A7 structure (the sc structure is a special 
case of the A7 structure), b) Intermediate structure. Compared to a) the unit cell of 
the sc lattice is elongated along the vector keeping the volume per atom constant. 
The new positions of the atoms are shown together with a cube of the sc lattice, 
c) The A7 structure of As at ambient pressure. Compared to b) the atomic planes 
perpendicular to a3 are displaced alternatingly in the a^ and — a3 directions due to a 
Peierls distortion. The cube still represents the sc structure of a), for reference. Below 
a) - c) projections of the atomic planes onto a^ are shown. 



2. Method 

To determine the phase of As under pressure before the laser excitation when the 
electrons are in their ground state we compared the enthalpies of the A7 and sc phases, 
which were computed in the following way. (i) We calculated total energies for a series 
of unit cell volumes with the computer program wien2k This is an all-electron 
full-potential linearized augmented plane wave code, which has been designed to make 
very accurate predictions for solids within the limitations of density functional theory 
[T7] . For the sc structure these calculations are straightforward, because there are no 
internal parameters. However, for the A7 structure the energy is a function of both 
volume and the internal parameters c/a and z described in the Introduction. In this 
case we minimized the total energy at each volume by optimizing the c/a ratio and 
the Peierls distortion parameter z. (ii) We fitted the total energy vs volume data for 
the sc phase to the Birch-Murnaghan equation of state [H] , which has been derived for 
systems with cubic symmetry. The fitting parameters of this equation of state are the 
minimal total energy Eq, the corresponding atomic volume Vq, the bulk modulus Bq, 
and the first derivative of the bulk modulus with respect to pressure B'q. (iii) We fitted 
the computed energy differences between the A7 and sc structures to 



We found that this form described our A7 data much better than when the total energies 
of the A7 structure are fitted directly to the Birch-Murnaghan equation of state [18]. 
(iv) We calculated the enthalpies H = E + PV for the sc and A7 structures from the 
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analytical fits for the total energy and -P = ~ (ff ) 



Laser-induced solid-solid phase transition in As under pressure 



4 



To describe As after laser excitation we use the following physical picture: The laser 
pulse creates electrons and holes, which undergo dephasing and collisions on a timescale 
that is much shorter than the typical time of ionic motion (~ 130 fs, based on the Aig 
optical F-point phonon frequency of As in its equilibrium structure). Therefore one can 
for all practical purposes assume that the excited carriers thermalize instantaneously. 
The effect of the excitation by an ultrashort laser pulse was thus simulated by heating 
the electrons. In this case, the phase stability is governed by the free energy F = E—TS, 
where T is the electronic temperature and S is the electronic entropy (To be sure, the 
free energy was also used for the ground state calculations, but since the electronic 
temperature was very low in that case and was used for smearing purposes only, we 
have ignored this point for the sake of clarity of presentation in the discussion above). 
Note that the entropic contribution of the ions to the free energy was neglected. This is 
because we assumed that both the laser pulse duration and the timescale on which the 
phase transition takes place are much shorter than the electron-ion interaction time, 
so that no substantial heating of the ions due to the laser occurs during the time of 
interest. Another consequence of the short timescales that we are looking at is that the 
system does not have time to expand, so that structural changes should be assumed 
to take place at constant volume. Therefore, we compared the free energies of the A7 
and sc structures to predict under which conditions a phase transition takes place in As 
under pressure. 

Details of our wien2k calculations were as follows. We used the generalized 
gradient approximation [19] for the exchange and correlation energy. Our basis consisted 
of plane waves with energies less than or equal to 20.25 Ry. Inside the so-called muffin- 
tin spheres around each atom the plane waves were augmented by atomic orbitals with a 
linearized energy dependence. To get a better convergence with respect to the number 
of plane waves included in the basis the 3d, 4s, and 4P states were described with a 
combination [20] of energy independent augmented plane waves and local orbitals. Also, 
additional local orbitals [21] were added for the 3d and the 4s states. The first Brillouin 
Zone of the A7 structure was sampled using a grid of 32 x 32 x 32 k points excluding the 
r point, which corresponds to 2992 irreducible k points. To rule out errors due to, for 
example, a different k-space sampling and to allow for a detailed comparison between 
the A7 and the sc structures, the latter structure was also computed as an A7 structure 
with the parameters c/a and z fixed to their sc values -\/6 and 0.25, respectively. 

The electronic ground state calculations were performed using temperature 
smearing (T = 1 mRy) to determine the electronic occupation numbers. For the 
laser-excited state an electronic temperature of T = 19 mRy was used. This choice 
corresponds to an absorbed laser energy of Eiaser ~ 2.8 mRy/atom or, equivalently, the 
creation of ~ 0.06 electron-hole pairs per unit cell (when the atomic volume V = 114.0 
Og/atom). Since we assumed an instantaneous thermalization of the excited carriers, 
our predicted results do not depend directly on the energy of the phonons, but only 
on the total absorbed energy per atom. The incident laser fiuence that leads, at the 
sample surface, to the above-mentioned value of Eiaser is, however, a function of both 
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Table 1. Fitting parameters for the electronic ground state (T = 1 niRy) and the 
laser-excited state (T = 19 mRy). Vq, Bq, and B'^ give the best fit of our sc data to 
the Birch-Murnaghan equation of state [T5]. T^, /3c, and A^. represent our best fit of 
the computed energy differences between the A7 and the sc structures to (IT]). 



T = 1 mRy T = 19 mRy 



Vo (ap/atom) 


137.479 


138.265 


(GPa) 


78.107 


76.241 


B', 


4.317 


4.341 


Vc (oo/atom) 


111.31 


115.54 




1.429 


1.611 


Ac (mRy/ atom) 


-48.63 


-59.37 



the reflectivity and the penetration depth of the laser light, which are dependent on the 
laser wavelength. A complication is that these dependencies might change due to an 
applied pressure and under the influence of the laser irradiation. Ignoring these latter 
points and using a laser wavelength of 800 nm we estimated our incident fluence to be 
~ 1.8 mJ/cm^ (for this wavelength the reflectivity is ~ 55% [22] and the penetration 
depth is ~ 23 nm, based on the experimentally determined complex dielectric constant 
of 16 + 2H [23]). Here we would further like to mention that we assumed the laser pulse 
duration to be much shorter than the timescale of ionic motion (~ 130 fs). 

3. Results 

Before turning to the energy vs volume curves, a few words about the equilibrium 
structure of Arsenic are at place, because this allows us to get a rough idea about the 
quality of the predictions that the GGA [19j we used is able to make. According to 
our calculations the ground state parameters of As in the A7 structure are c/a = 2.86, 
z = 0.226, and V = 154 Op/atom. This is in reasonable agreement with the measured 
values c/a = 2.78, z = 0.228, and V = 144 a|]/atom jl5]. As was to be expected, 
our atomic volume V is larger than the experimental value [15] by about 7%, which 
is consistent with the general trend of the GGA to overestimate bond lengths, and in 
agreement with the value of 154 Og/atom recently obtained from GGA calculations by 
Shang and co-workers [Ij. From calculations for smaller atomic volumes we concluded 
that our too large values for the magnitude of the Peierls distortion (0.25 — z) and for 
c/a can in part but not entirely be understood as a consequence of the overestimated 
unit cell volume. 

The total energies of the sc and A7 structures as a function of atomic volume are 
shown in figure |2} These data are for the electronic ground state, i.e., before laser 
excitation. The parameters obtained from the fits to the Birch-Murnaghan equation of 
state and ([T]), which we made using the data points plotted in figure |2] and the inset 
of figure [2] respectively, are given in table [T] From figure |2] it is apparent that our fits 
follow the computed data very closely. Indeed, the root mean square of the residuals 




Figure 2. Total ground state energy vs volume for sc As and As in the A7 structure. 
The solid curve shows the Birch- Murnaghan fit to the data for the sc structure. The 
inset shows the energy difference AE — Eat — Egc (plusses) and the fit to ([T]) (solid 
curve). The volume Vc at which the A7 structure becomes energetically favored over 
the sc one is indicated. 




Pressure (GPa) 



Figure 3. Enthalpy difference AH = H^j — Hgc as a function of pressure. The inset 
shows the indicated region enlarged. The transition from the A7 to the sc structure 
takes place at a pressure of 26.3 GPa. 



was only 0.003 mRy/atom for the Birch-Murnaghan fit and 0.005 niRy/atom for the fit 
to 0. 

The enthalpy of the A7 structure relative to the enthalpy of the sc structure as a 
function of pressure is plotted in figure |3} In a narrow pressure range multiple solutions 
exist for the enthalpy of the A7 structure, as is shown in the inset of figure |3| Of course, 
the solution with the lowest enthalpy is the stable one. A consequence of these multiple 
solutions is that the enthalpies of the A7 and sc structures cross and that the transition 
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Figure 4. a) the Peierls parameter z and b) the ratio c/a as a function of volume 
before (T = 1 mRy) and after laser excitation (T = 19 mRy). The inset of b) shows a 
selected region enlarged. Arrows in a) and the inset of b) indicate the possibility of a 
laser-induced phase transition at constant volume from the A7 to the sc structure. 



point from the A7 to the sc structure can thus easily be determined. From figure |3] we see 
that the transition takes place at 26.3 GPa, in reasonable agreement with earlier GGA 
calculations by Haussermann and co-workers |21] , who found the transition pressure to 
be 28 GPa. Another consequence of the multiple solutions seen in the inset of figure 
[3] is that the unstable solutions, corresponding to a narrow range of possible atomic 
volumes, are not realized for any given pressure, which means that there is a small 
volume change at the transition, which according to our calculations is AV/Vc = 0.4%. 
Both the transition pressure and the volume change are in nearly perfect agreement with 
the experimental values of Beister and co-workers p] , who have found the transition at 
a pressure of 25 ± 1 GPa accompanied by a volume change AV/Vc ^ 0.6%. We did not 
find a discontinuity in c/a at the transition, which according to Beister and co-workers 
[2] jumps from 2.49 to ~ 2.45. However, according to our calculations z changes 
discontinuously from 0.246 to 0.25. 

We now describe the effect of an ultrashort laser pulse. As we have already argued 
in section [2]the timescale that we are interested in is too short for the system to expand. 
Therefore we present results at constant volume, and we only briefiy discuss the effect of 
pressure at the end of this paragraph. Our optimized c/a ratio and Peierls parameter z 
as a function of volume for As in the A7 structure are shown in figure |4j It can be seen 
that the pressure-induced phase transition from the A7 to the sc structure occurs at a 
larger volume when the electrons have been heated with a laser compared to when they 
are in their ground state. As a consequence, it is possible to induce a phase transition at 
constant volume from the A7 to the sc structure in As under pressure. We have indicated 
this possibility for a specific volume (V^ = 114.0 Og/atom) by two arrows in figure |4| 
The corresponding potential energy surfaces before and after laser excitation are shown 
in figure |5| When the electrons are in the ground state this volume corresponds to 
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Figure 5. Potential energy surface contour plot of As a) before and b) after laser- 
excitation as a function of the internal parameters z and c/a at the volume V = 114.0 
ap/atom. Whereas the energy minimum of As before laser-excitation is relatively sharp 
(the two minima shown are equivalent), the minimum after laser-excitation is broader. 
The minimum in b) is located at the sc values of z and c/a (0.25 and sa 2.45, 
respectively) . 

an applied calculated pressure of 23.8 GPa. The optimized ratio c/a = 2.48 and the 
Peierls parameter z = 0.242 indicate that the system is in the A7 phase. After the 
laser excitation the optimal c/a and z become 2.45 and 0.25, respectively, signifying 
that As has made the transition to the sc structure. At this given volume, the internal 
pressure in laser-excited As has dropped to 22.4 GPa. If the system were given enough 
time to adjust to the apphed pressure of 23.8 GPa, the volume would, according to our 
calculations, decrease to 113.1 ajj/atom, in which case the simple cubic phase would 
clearly still be the most stable one (see figure |4]). 

Comparing the free energies of the A7 and sc structures it follows that the A7 to sc 
transition in As can, in principle, be induced for atomic volumes up to < = 115.54 
ag/atom (T = 19 mRy data in table [T]), which corresponds to an applied pressure of 
> 22.1 GPa. However, near the transition point the potential energy surface of As 
becomes very flat (see figure [s]), and as a consequence the uncertainty in the optimized 
c/a and z values becomes relatively large, which explains why in figure |4] the c/a ratio 
and z Peierls parameter start to deviate from their sc values at volumes less than 
Vc = 115.54 ttg/atom (but greater than V = 114.0 al/atom). So, in order to be on the 
safe side and to avoid any ambiguity about the final state after the laser excitation, the 
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lowest pressure for which we feel confident to predict an unambiguous transition to the 
sc structure (as defined by its c/a and z values) is the value of 23.8 GPa derived in the 
previous paragraph. 

As can be noted in figure |4](b), the predicted phase transition involves a change in 
c/a. Since the volume remains constant, in a polycrystalline sample there is no average 
effect and we expect that the change in c/a takes place on a subpicosecond timescale. 
In a monocrystal this change might be limited by the sound velocity. 

4. Conclusion 

Summarizing, we demonstrated that it is possible to induce a solid-solid phase transition 
in As under pressure using an ultrashort laser pulse. In particular, for an absorbed laser 
energy Eiaser ~ 2.8 mRy/atom, we predicted that the phase transition will be induced 
under an applied pressure as low as > 23.8 GPa. For higher absorbed energies we 
expect that the transition could be induced under an even lower applied pressure. This 
transition might be observed using time-resolved x-ray diffraction experiments. 
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